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Abstract 

We analyse the low lying spectrum of a model of excitons in carbon 
nanotubes. Consider two particles with opposite charges and a Coulomb 
self-interaction, placed on an infinitely long cylinder. If the cylinder radius 
becomes small, the low lying spectrum of their relative motion is well de- 
scribed by a one-dimensional effective Hamiltonian which is exactly solvable. 
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1 Introduction 



In order to understand the quantum mechanics of carbon nanotubes, one has to 
reconsider many classical problems in which the systems are restricted to low di- 
mensional configuration spaces. The effects induced by these special shapes are sig- 
nificant. For example, optical properties and electrical conductivity in nanotubes 
and nanowires are highly influenced by their geometry. 

In a periodic structure, bands of allowed and forbidden energies are char- 
acteristic for non-interacting electrons. When applying an external perturbation, 
such as light, electrons can only absorb the amount of energy which allows them 
to jump from an occupied level of energy to a free one. In the particular case of a 
semi-conductor, at low temperatures, the band of energies are either full (valence 
bands) or empty (conduction bands). So the electron must absorb a fairly large 
amount of energy to jump to the conduction band. 

When the self-interaction is also considered, the mathematical problem of the 
optical response becomes very difficult, and there are not many rigorous results 
in this direction. Here is a sketchy description of what physicists generally do (see 
for example the book of Fetter and Walecka |FWj : 

1. Work in the grand-canonical ensemble, at zero temperature, and the Fermi 
energy Ep is in the middle of an energy gap; 

2. Switch to an electron-hole representation, via a Bogoliubov unitary transfor- 
mation. The new non-interacting ground state is the tensor product of two 
vacua. If before an excited state meant that an electron was promoted from 
an occupied energy state from below Ep to an empty state above Ep, in the 
electron-hole representation it simply means that an electron-hole pair was 
created; 

3. Try to diagonalize in one way or the other the true, self-interacting many- 
body Hamiltonian by restricting it to a certain subspace of "physically rele- 
vant" excited states; this usually amounts to formulate a Hartree-Fock prob- 
lem; 

4. Try to obtain an effective one-body Hamiltonian, whose spectrum approxi- 
mates in some way the original problem in a neighborhood of Ep; 

5. The bound states of this effective one-body operator are called excitons. 
They describe virtual, not real particles; 

6. Use the exciton energies to correct the optical response formulas derived in 
the non-interacting case. 

Now this one-body effective Hamiltonian is a complicated object in general. If one 
makes a number of further assumptions like: 
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1. There is only one conduction band above Ep and only one valence band 
below Ep; 

2. The dispersion law of these two energy bands is replaced with paraboloids 
(effective mass approximation), 

then this one body effective Hamiltonian is nothing but the one describing the 
relative motion of a positively charged particle (a hole) and a negatively charged 
particle (electron), interacting through an eventually screened Coulomb potential. 
Let us stress that this procedure is generally accepted as physically sound in the 
case when the crystal is periodic in all three dimensions. 

If a special geometry is imposed (i.e. the electrons are confined on long and 
thin cylinders made out of carbon atoms) then the above procedure has to be 
completely reconsidered. The problem is even more complicated, because two di- 
mensions are on a torus and the band structure only arises from the longitudinal 
variable. Even the position of the Fermi level moves when the radius of the cylinder 
varies. 

It has been argued by physicists PedlL IPed2| that one can still write down 
a Hartree-Fock type eigenvalue problem which describes the excited states near 
the Fermi level. This operator is a two-body one, which does not in general allow 
a complete reduction of its mass center. A mathematically sound formulation of 
this Hartree-Fock problem would be of certain interest, but it is not what we do 
in this paper. We rather perform the spectral analysis of an operator which has 
been conjectured by physicists as being the relevant one. 

The main point in investigating these low dimensional structures, is that the 
interaction between electrons is enhanced and gives much stronger exciton effects 
than in three dimensions. This means that some new energy states appear deep 
inside the forbidden energy band. The smaller the radius, the more important 
these new energies are. That is because they allow photons with much smaller 
energy than the band gap to be absorbed into the material. 

We will therefore consider two charged particles living on a cylinder and 
interacting through an attractive Coulomb potential. As we have already pointed 
out above, this operator models an effective Hamiltonian for excitons in carbon 
nanotubes, according to |Pedl| and |Ped2j . Remember that it has nothing to do with 
real particles living on the nanotube, the exciton being just a mathematical artifact 
describing virtual particles. We hope that our current results could also describe 
excitons living in nanowires |Aki| . or two particles in a strong magnetic field as 
in |BU| . Let us mention that our paper is an improvement and a continuation of 
a previous work done in |CDP| . 

2 The mathematical model 

Our configuration space is a cylinder of radius r and infinite length, space denoted 
by Cr = M X rS^ , being the unit circle. The coordinates on the cylinder are 
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{x,y) e (R X rS^) where x is the variable along the tube axis and y is the transverse 
coordinate. 

The two virtual particles live in the Hilbert space L^(Cr x C^). We formally 
consider the Hamiltonian 

H^ = -h'(^ + p^ + p^ + p^)- V^{x^ - X2, yi - y,), (2.1) 
\2mi 2m2 2mi 2m2 / 

where 

Vr{x, y) := = (2.2) 

eJx'^ +4r2 sin^ (^) 

{xi,yi) axe the coordinates on the cylinder of the two charged particles, rrij their 
masses, and e, their charges which are of opposite charge so that Vr > 0. Here e is 
the electric permittivity of the material. In the sequel we will set H = e = 1. The 
potential y is the three dimensional Coulomb potential simply restricted to the 
cylinder. We justify the expression of V" by Pythagora's theorem. The cylinder is 
embedded in R^. The distance p from one particle to the other in M.^ is: 

p = [xi - X2) + 4r^ sm ( — — — 

where |2r sin ( ^^2r'^ ) I is the length of the chord joining two points of coordinate 

yi and y2 on the circle. 

Now consider the space 

Vo = {f eC^iCrXCr): Vai,2,/3i,2,7i,2eN (2.3) 

\x^^xrD^lD^,lD^^lD^^lf{xuyi,X2,y2)\<C„,0,j } 

of " Schwartz functions along x" and smooth and 27rr-periodic along y. Clearly Dq 
is dense in the Sobolev space H^{C^). Let us define on Vq the kinetic quadratic 
form 

torn = ^ iWd.Ml' + WdyM') + ^ iWd.M' + WdyM') (2.4) 
and the quadratic form associated with the Coulomb potential on the cylinder: 

tvm = \\VVri^f, (2.5) 

/ Vr{xi-X2,yi-y2)\ip{xi,X2,yi,y2)\'^ dxidx2dyidx2<oo ]. 
Finally, define the form 

tn '■= to — W on the domain Dq. (2.6) 
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The scsquilincar form induced by to is densely defined, closable, symmetric, non- 
negative, and its corresponding selfadjoint operator Hq is —-^^Ai — with 
periodic boundary conditions in the y variables. Its form domain is H^{C^), and 
is essentially self-adjoint on Vq. 



2.1 Center of mass separation in the longitudinal direction 

We introduce the total mass M := nii + m2 and the reduced mass /x := ^^^2 " 
Denote as usual with 

( Y _ mixi+m2X2 ™ _ _ 

J ^ ~ mi+m2 ; ^ - ^2, 

\ Y = y2, y = yi - y2 

then, 

X,=X + ^X, X2=X- ^x, 



yi=y + Y, y2 = Y 

Unfortunately, for the yi, y2 variables we cannot use Jacobi coordinates because 
the transformation does not leave invariant the domain of the Laplacian (the pe- 
riodic boundary conditions are not preserved). That is why we use atomic coordi- 
nates y and Y instead. In these new coordinates, the total Hilbert space splits in 
a tensor product L'^{M.) L'^[M.x (rS^)"^]. More precisely, if we denote by 

U-^ : L''iR)<g)L^[R x {rS^)^] ^ L^Cl), [U'^ !]{xi,yi,X2,y2) = .f(X,x,Y,y) 

then it is quite standard to show that after this variable change we can separate 
away the X variable and for f,gG UIVq] we get: 

tH{U-'f,U-'g) = ^{dxf,dx9) + ^^{9Yf,dYg) (2.7) 
1 



{9yf,dYg) - {^/Vr{x,y)f, ^Vr{x,y)g). 

Note that the subset J7[I?o] has the same properties as Vq but in the new variables. 
Therefore we can concentrate on the reduced form 

Hf,g) = ^^{dYf,dY9) (2.8) 
^ 2^^^^-^' ^^^^ ^ ^^^^ 



- —{9yf,dYg) - WVr{x,y)f, ^/Vr{x,y)g), 



m2 
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densely defined on smooth enough functions in L'^[R x (rS'^)^], dc }caying along 
the X variable. Consider the decomposition L^[]R x (rS^)"^] = ^,^^j^L^{R x rS^) 
implemented by the Fourier series 



iYi 



fix,y,Y) = V -j==fk,r{x,y) 



where 



1 rZTTr 

fk,r{x, y) = -= / f{x, y, y)e-'^*dy 
V27rr Jo 

Then for our form we get: 



th=:^th„ (2.9) 



feez 



where t^^. is 



th,{f,9) = ^{d.f,d,g) + ^{dyf,dyg) + ^^(Lg) (2.10) 
ik 



defined on the domain 

£k - {/ € C°°{m X rT) \x''D^:Dlf(x,y)\ < ^/j^, Va,/3,7 € N} . 

Now remember that we are only interested in the low lying spectrum of our original 
operator. We will now show that for small r, only thg contributes to the bottom 
of the spectrum. Indeed, let us concentrate on the operator 

2iJb " 2m2r^ miv^ 

defined on rS^ with periodic boundary conditions. Via the discrete Fourier trans- 
formation it is unitarily equivalent to: 

(2 111 



where 



- \^.^}^(K -¥^](l)- (2.12) 



2/ir2 2m2r2 m2r'^ 2r2 V";^ 7^ / 
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A simple calculation shows that both eigenvalues of the above matrix are positive; 
denote with A_ the smaller one. Then the operator in (|2.11() obeys: 



mf±P! + _^! ^^>_^a2 A_fc^ ,2 131 

^ I 2^ r2 ^ 2m2r2 msr^ / " 2 ^ 2r2 • ^ ' 



Using this in H2.10|l we obtain the inequality 



tk, if, f) > min {1, flX^} ^mf, dj) + {dyf, dyf)] 

-iHAfJ) + ^\\f\? (2.14) 

where thg is obviously defined by the previous line. Now one of the results obtained 
in this paper will be that the spectrum of the self-adjoint operator associated to 
a form like th„ is bounded from below by a numerical constant times — (ln(r))2. 
Hence if fc ^ and r is small enough, all th^. will be positive and only thg will 
contribute to the negative part of the spectrum. 



2.2 The self-adjointness problem 

Due to H2.7|l . H2.8|l . (|2.9|l and H2.14|l . it is clear that it is enough to concentrate 
on tho ■ If we can prove that it is bounded from below, then all other forms with 
fc ^ will also have this property, and the total Hamiltonian will be a direct sum 
of Friederichs' extensions. Because we can anyway scale the masses and charges 
away, and in order to simplify the notation, let us consider the sesquilinear form: 

t^{f,9) l[{d.f,d^g) + {dyf,dyg)] - VM^)9) 

^■.W,g)~W{f,g) (2.15) 

on the domain 

£ = {fe C^{^ X rS^) \x'-DPDlf{x,y)\ < C^p, a,/3,7 € N} , 

now where Vr is as in (|2.2|l but with 6162 = —1. 

We will now construct a self-adjoint operator out of this form. 

Proposition 2.1. The operator Vr is relatively compact in the form sense with 
respect to the operator —^{d^ + dy) with form domain 7i^(Cr). Thus the form tn 

defines a self adjoint operator whose form domain is Ti}(Cr), and a^ssiH'') = 
[0,00). 
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Proof. Wc identify the cylinder Cr with the strip R x [— ttt, ttt]. For every a > 
we define 71°' {Cr) to be the set of all functions which (at least formally) can be 
expressed as: 

/(x,2/) = I e^^^+™«/Vn.(?') dp, 

J2 [{1 + bl'" + \mfn\Lip)\'dp < oo. (2.16) 

Let X be the characteristic function of the interval [— r/2,r/2]. Since near 
the boundary of the strip Vr{x, y) • (1 — xiu)) is bounded, we only have to look at 
V{x, y) := V{x, y) ■ xiu)- Then we can find a constant C such that everywhere in 
Cr we have 

\V{x,y)\< S^ . 

Denote with p := \J x?- + . Choose a function Xi G C^(]R) with support in 
(— 3r/2, 3r/2), such that xXi = X- Then the operator of multiplication by xi is 
bounded from H"{Cr) to H"(M^) and vice versa, because it does not touch the 
boundary (the proof of this fact is standard). Moreover, if —A is the operator 
associated to to, then we have 

H^iCr) = (-A + ly/^L'^iCr). 

Note that £ is dense in any TL°'{Cr)- Moreover, for every G £ we have 

KV', Vi>)L2^Cr)\ < C{xii>, ^XiV')l^(r=)- (2.17) 
We have that XiV' € »5(M^). Then we can write 

(Xi^, (l/p)XiV'>L^(M^) = / / Xi(psin(0))V^(p, 0) . ^jj{p, 9)dpde 
Jo Jo 

and after integration by parts in the radial integral we obtain 



(Xi^, (1/p)xiV')l2(i:2) 

/-•27V nOO 

= -/ / {dp[xiipMO))^ip,o)-xi{psm{om{p,e)]}pdpde. (2.18) 

Ja Jo 

Then using the estimate \dp{xi'fp)\ < |V(xiV')|! and with the Cauchy-Schwarz 
inequality: 

(XlV", (1/P)X1V'>L2(R2) < 2||xiV'IIl2(R2)||V(X1'0)||l2(r2) 

<const||V||L2(c,)||V||7ii(c.). (2.19) 
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Now for an arbitrarily small e > we have 

l(V',t^^>L^(c.)l < (Ci/e)||V^||i.(c.)+^iei~o(V',V'), 

where Ci is just a numerical constant. The density of E in T-L^(Cr) finishes the 
proof of relative boundedness, and we can define as the Frederichs extension. 

Until now we have shown in an elementary way that A + is 

bounded, but one can do much better than that. In jEo| it has been proved a two 
dimensional version of an inequality of Kato, which states the following: 

rfi/4)* , 

(V-, |x|-V)l^(r^) < ^^^2 \/^V')l^(r^)- (2.20) 

This inequality immediately implies that ^/Vr '■ Ti^^^{Cr) L'^{Cr) is bounded. 
Now let us show that the operator -\/T^^^^^(Cr) L^(Cr) is compact. We will 
in fact prove the sufficient condition that the operator T := |x|~^/^(— A + 
defined on L^(R^) is compact. 

Indeed, let us denote by Xn the characteristic function of the ball of radius 
n > 0, centered at the origin in R^. Then we can write: 

T = Xn(x)Tx„(-A) + [(1 - Xn)(x)]r + Xn(x)r[(l - x„)(-A)]. 

First, the operator Xn(x)Txn(— A) is Hilbert-Schmidt (its integral kernel is an 
L^(]R'*) function), thus compact. Second, the sequence of operators [(1 — Xn)(x)]r 
converges in norm to zero. Third, the sequence Xra(x)T[(l — Xn){—^)] can be 
expressed in the following way: 

X„(x)r[l - x„(-A)] = {x„(x)|xri/^(-A + l)-i/4}(_A + l)-i/4[(i _ Xn)(-A)], 

where the first factor is uniformly bounded in n, while the second one converges 
in norm to zero. Thus T can be approximated in operator norm with a sequence 
of compact operators, hence it is compact. 

Therefore 1^ is a relatively compact form perturbation to —A, hence the 
essential spectrum is stable, and the proof is over. □ 

2.3 An effective operator for the low lying spectrum 

We will show in this section that at small r, the negative spectrum of can be 
determined by studying a one dimensional effective operator -ff^ff . It is natural to 
expect that the high transverse modes do not contribute much to the low region 
of the spectrum. 

First, we separate into differejot parts taking advantage of the cylindrical 
geometry, that is to say, we represent as a sum of orthogonal transverse modes 
using the periodic boundary conditions along the circumference of the cylinder. 
Second, we analyze which part is relevant when the radius tends to zero. 
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We recall that is formally given by W' = — ^ — Vr in the space 
L^(Cr) ~ L^(IR) (g) L^(r5'i). The domain contains all ^ G 'H^{Cr) with the property 
that in distribution sense we have 

A. 



2 2 Vrj^eL^Cr). (2.21) 

This does not mean that the domain is H^{Cr) because Vr is too singular at the 
origin. 

Our problem has two degrees of freedom. We consider the orthonormal basis 
of eigenvectors of — with domain ■Wpg^((— 7rr, 7rr)) ~ 7{^{rS^). Here, the Sobolev 
space Wpgr((— Trr, 7rr)) denotes functions which are 27rr-periodic with first and 
second derivatives in the distribution sense in L^. We can write 

= £ Kiii 

n— — oo 

where the one dimensional projectors are defined by 

n; = (x;, •)x;, xl{y) = ;7^e-^ and = n e z. 

We now introduce a family of orthogonal projectors 

n:;:=i®n;, (2.22) 

which project from L^(Cr) into what we call the n*'* transverse mode. The operator 
can be split as follows: 

Hr = Y,tKJI^tKn =: ® (^^n' '^^n) + E ^n,m ® ((xL, ■)xl), (2-23) 

where the sum is a direct sum, since the projectors are orthogonal. By a natural 
unitary identification, we can work in a new Hilbert space: 

H = l'[Z-h\^)], H 9 V = {V'ninez, e L2(R). (2.24) 

Therefore our original operator is an infinite matrix now, {Hn,m}n,mei. whose 
elements are operators in L^(R). 

lin^m, the only contribution comes from Vr, and the corresponding oper- 
ator is a multiplication operator given by (x ^ 0): 

27rr J_^r 

If n = m, then the corresponding diagonal element is given by the operator: 



10 



where V^j is deduced from when m = n and is given by 

V:s{x) = ^ £^ V{x,y)dy. (2.27) 

Finally, let us introduce a last notation for what will be our effective one-dimensional 
comparison operator: 

^eff:=-^i-Kff(2^) (2.28) 



2dx2 

and note that 



n2 



Hl = Hl^ + ^. (2.29) 

One can see that for n ^ 0, the diagonal entries of our infinite operator valued 
matrix are pushed up by a term proportional with Thus a natural candidate 
for a comparison operator for the negative spectrum of H'^ is H^q- In the next 
section we will perform a careful study of this operator. 

3 Spectral analysis of 

We now want to study the spectrum of the operator H^^ when r becomes small. 
We recall that: 

^- = -^^-^-(^) 

where 

1 r'' 1 

KffW = IT- / ; dy 
27rr Ja;2+4r2sin2^ 



with form domain Q{H^ff) = H^(M). We are going to use perturbation theory 
around r = 0, which will turn out to be quite a singular limit. The strategy is to 
approximate the form associated to the potential VJg{x) around r = by another 
quadratic form which provides a solvable approximation. 

Let us define the sesquilinear form on S{R) (later on we will show that it is 
bounded onH^' 



poo pO 

Coif, 9)-=- lii(2x) • [f{x)g{x)]'dx + / ln(-2x) • [f{x)g{x)]' dx 

Jo J-oo 

= - / ln(a;) • [f^g{x)]'dx + / ln(-a;) • [f{^g{x)]'dx + ln(4)7(0)5(0) 

Jo J-oo 

= |fp^ + ln(4)5|(75). (3.1) 

The symbol fp means the finite part in the sense of Hadamard, while 5 is the Dirac 
distribution. Note that up to an integration by parts, and for functions supported 
away from zero, we have Co{f,g) = (/, j^g). 
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The main result of this subsection is contained in the following proposition: 
Proposition 3.1. For r < 1 and for every f G W^(M) we have the estimate 



(/,Kff/> = -21n(r) 1/(0)1^ + Co(/, /)+ O (r^) 



= -21n(r/2) |/(0)|^ + 



\x\ 



llwi(R)- 



(3.2) 



Proof. The argument is a bit long, and we split it in several lemmas. Let us start 
by listing some of the properties of V^f^. First note that it scales like a "delta 
function", i.e. it is homogeneous of order —1: 



(3.3) 



The next observation is that due to the integral with respect to y it is much 
less singular than Vri 

Lemma 3.2. The behavior ofV^g{x) is logarithmic at 0. 

Proof. There exists a constant c > large enough such that for every |?/| < f and 
a; e K we have 

1 x"^ 



Thus we can integrate and obtain V^g(a;^ 



■In{\x\) + 0{1). 



□ 



We now define on K a comparison function Yr{x) := -^7=5==^. and we also 

denote by Yr the associated quadratic form defined on ^.^{M.). For the following, 
let us recall the classical Sobolev estimate in one dimension: 

||/IU<^(||/'|| + 11/11) = ^||/lk,(M). (3.4) 

Lemma 3.3. We have the following properties: 

(i) Kff > Yi, V^six) - Yi{x) = 0(|x|-5) for \x\ > 10, and 

||Kff-illUHM)=M4); (3.5) 

(ii) (/, (y^ - Yr)f) = ln(4) |/(0)p + 0(rt)||/||2^,(«, Vr < 1. 
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Proof, (i). To show that V^g > Yi, one uses | sin(-)| < 1. The second estimate for 
Ixl > 10 follows from: 



V,],{x)-Y,{x)^ — 



4sin2f 



dy 



\x\ 2tt 



1 r 2 sin^ ^ 



dy + 0{x-^)~l + — +0{x-^) 



\x\ > 10. 



(3.6) 



Before computing the norm of (|3.5|l . let us notice that none of the terms is in 
L^. We first integrate with respect to a:, and then over y, and get: 
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K^ff - r^lUnR) = / ln(siny)d2/ - 21n(2) 



(3.7) 



thus 13.5|l is proved. 

Let us now prove (ii). We have due to the scaling properties: 

(/, (K'ff - Yr)f) (K'ff - Y,){x/r)\f{x)\'dx 



= / {V^s~Y^){x)\f{rx)\'dx. 



(3.8) 

Then, we subtract the term \\V^g - Yi\\li ■ |/(0)P ln(4) • |/(0)|2 which gives 
(/,(K'ff-i;)/)-hi(4).|/(0)p = 

{V^s - Y^){x) {[f{rx) - /(0)]7M + [7M - 7(0)]/(0)| dx. (3.9) 

Let a G (0, 1) a real number. We split the above integral in two regions: < r^" 
and > r^". We have, using (|3.4() : 



/ (K'ff - Y,){x)[f{rx) - fmf{rx)dx 

J|£c|>r-° 

<2||/||?,i / \{V^g-Y,){x)\dx 

J|a;|>r-° 

< 11/11?.^ / i'^^ ifr-">10 

J\x\>r-'^ l-^l 



(3.10) 



„4q\ IIJ-II2 



<0(r-"). 11/11^.. 



(3.11) 
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For the region < a; < r " (and similarly for the other one), we can write: 

/ [V^s - Yi){x)\f{rx) - /(0)| \f{rx)\dx (3.12) 

J0<x<r-° 



f{t)dt 



\f{rx)\dx 







< / (Kff->^i)(^) 
<llK'ff-Fi|U..||/IU~. / \f\t)\dt 

Jo 

and the Cauchy-Schwarz inequality yields: 

|/'(i)Mi<r^||/||„i. 
Then we set a as the solution of (1 — a)/2 = Aa which gives a = h. □ 



We now concentrate ourselves on Yr when r is small. For the next two lemmas, 
we need to introduce the the following characteristic function: 

= { I othelwist ■ (^-^^^ 
Then we have the following lemma: 

Lemma 3.4. Consider the self-adjoint operator of multiplication by ln(| • |)x 
defined on its natural domain in L^(R). This operator is relatively bounded to 
Px '■= —id/dx, with relative bound zero. 

Proof. Indeed, ln(|a;|)x(|2:|) {px + iA)^^, A > 1 is Hilbcrt-Schmidt since we have, 
from XI.3]: 



II ln(| • |)x (p. + iX)-'\\hs < const • || ln(| • |)x||lHI(- + «A)-i|| 



L2 



const 

< ,3.14, 

Note that by a similar argument as the one in (|3.4|) we get the estimate: 



,N 1,1 const 
||(p, + zA)-1||l2_L- ^ (3-15) 



Then a standard argument finishes the proof. □ 
We can now characterize the form Cq introduced in 1)3. l|l : 
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Lemma 3.5. The quadratic form induced by Co admits a continuous extension 
to 7i^(M). Moreover, Cq is infinitesimally form bounded with respect to the form 
associated to = —d^/dx^. 

Proof. Fix some e G (0, 1). Then for every / G iS(M) we can write: 

Coif J) = - / ln(2:E) • / ln(-2:E) • {dM\'')ix)dx 



\n{-2x) ■ idMnix)dx - I H2x) ■ {d,\f\')ix)dx 

1 



ln(2e) + + / _ ^ ^ • Ifix^dx. (3.16) 



First we have 



/ ^■\f{x)\'dx<hfr. 

jR\[-e,e] Fl ^ 

Then using (|3.15|) we have 

snp\f{t)\<^\\ip, + ^X)f\l (3.17) 
tea V A 

which takes care of the terms containing /(±e). 

The remaining two integrals can be treated with the following argument: 



<2 / |ln(2x)|-|/(a;)|-|/'(a;)|da: 



\n{2x) ■ {dM\^){x)dx 

< 2||x ln(2 I . |)/|| ||/||„. < ^\\ip, + ^X)f\\ (3.18) 

where in the second inequality we used the Cauchy inequality, while in the third 
inequality we used (|3.14(l . 

These estimates allow us to find two constants A and B (growing when A 
grows, such that: 

|Co(/,/)|<^||/'|P+A||/||||/'||+i?||/|r 

But we use the inequality ||/|| ||/'|| < ^\\f'\\^ + XA\\f\\^ , which finally allows us 
to say that for any < a < 1 we can find & > such that 

|Co(/,/)| <a|b./|p + 5||/|p, (3.19) 

and the proof is over. □ 
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The final ingredient in proving Proposition 13.11 is contained in the following 
lemma: 

Lemma 3.6. Recall that Yr{x) — {x^ + 4r^)^^/^. Then for every r < 1, and for 
every f € S{M.), we have the estimate: 

if, Yrf) = -2 ln(2r) • |/(0)|2 + Co{f, f) + (r^) • WffnHm- 
Proof. Integrating by parts we obtain: 
(/,y,/) =-21n(2r)|/(0)|2 



ln(a; + ^/x^+Ar^) ■ [f'{x)f{x) + f{x)f'{x)]dx 



and: 



ln{-x + \/x2 +4r2) • [f{x)f{x) + f{x)f'{x)]dx 
(/,>;/> -Co(/,/) + 21n(2r) -1/(0)12 = 



(3.20) 



(3.21) 



[\n{x + y/x^+4r^) - \n{2x)] {d^\f{x)\^)dx 



+ / [ln(-a; + V^J^+iH) - ln(~2a;)] {d:r\f{x)\^)dx. 



The idea is to show that the last two integrals are small when r is small. We 
only consider the first integral, since the argument is completely analogous for the 
second one. 

Fix some < a < 1 (its optimal value will be chosen later) , and assume that 
r is small enough such that r^^" < 1/10. We split the domain of integration into 
two regions: one in which a; > r", and the other one where x < r". For the first 
region we have: 



ln(a; + \/aF+ir^) - ln{2x) = In ■ 



1 



1 



= 0(rVa;2). 



(3.22) 



Then by integration, and using H3.4|l together with the Cauchy inequality, we get: 



[ln(.T + +4r2) - \n{2x)]f{x)f{x)dx 



< 



< 



[ln{x + Va;2 +4r2) - ln(2x)]/'(a;)da; 



[ln(a; + Vx^ + 4r2) - \n{2x)]^dx \\f\\l^^ 



<0(r^-^")||/||^i. 



(3.23) 
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For the region where x < r", we use the monotonicity of the logarithm and write: 



I ln(a; + V^^+I^) - ln(2a;)| < | ln(a; + y/x^+Ar^)] + \ ln(2x)| 

< |ln(2r)| + |ln(2a;)|. (3.24) 

Then we can write 



[ln(a; + y/x^+Ar^) - \n{2x)]f{x)f{x)dx 



< const ■ \\f\\^^ [ln{2x)+ln{2r)]^dxj = 0{r^ \lnr\) ■ \\f\\^^. (3.25) 

Comparing (|3.23|) and (|3.25|l . we see that we can take a arbitrarily close to 1. In 
particular, we can find some a such that 2 — 3a/2 > 4/9 and a/2 > 4/9 and we 
are done. □ 

We can now conclude the proof of Proposition 13.11 by putting together the 
estimates from Lemma 13.31 and Lemma 13.61 

□ 



Remark. One can improve the exponent 4/9 in the error estimate (|3.2|l . and 
obtain 1/2 instead. One observes that the Fourier transform V^g(p) can be exactly 
computed in terms of modified Bessel functions, and then one expands it near p = 
identifying the Fourier transforms of fpj^ and S. Then the error's Hilbert-Schmidt 

norm is estimated and shown to be of order r^^^. A related problem treated with 
this method can be found in [5^2] . 

3.1 A solvable comparison operator He 

Remember that we are interested in the negative spectrum of i?cff' operator given 
in H2.28|l . Lemmas 13 . 21 and l3 . 41 1 ell us that its operator domain is 7i'^(R), while the 
form domain is 7i^(]R). Proposition 13. II see (|3.2|l . indicates that a good approxi- 
mation for H^g at small r would be the operator formally defined as 

Hc:=lpl + 2\n{r/2)S--^^. (3.26) 

Of course, as it is written above He makes no sense. The correct definition of 
He can be found in the Appendix A of |BD2| in a more general setting. For the 
comfort of the reader we give below a version of this definition adapted to our 
simpler situation. He has to be understood in the following way: consider the 
sesquilinear form on S{R) given by 

tc{f,g) ■■= l{f,g') + 21n(r) 7(0)g(0) - Co(/,ff). (3.27) 
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A standard consequence of H3.17|l and (|3.19|l is that the quadratic form associated 
with tc is closable, bounded from below, and the domain of its closure is H^{M.). 
Then He is the self-adjoint operator generated by tc, and its operator domain 
Dc is characterized by: 

Dc ■■^{ipeHH^) ■■ \tc{(t>,i^)\< const \\(t>\\, V<?i G 7i^(M)}. (3.28) 

Moreover, if ?/; e Dc, then wc have the equality: 

tc{^,^) = {(t>,Hc^), y^enHm- (3.29) 

Another representation for i/j G Dc is that there exists = Hcip G i^(M) such 
that the distribution ip" is a regular distribution on R \ {0} and is given by: 

i^"{x)^~2^^,{x)~2U{x). (3.30) 
\x\ 

One important consequence is that tp' E Ti.^{M. \ {0}), and ijj' is continuous on 
M\{0}. 

Let us now introduce the parity operators P+ and P_ 

P±:H'iR)^H\R), f{x)^{P±f)ix)^l^^l^J^. (3.31) 

We have that P+ + P- = 1. We will call RanP+ the even sector and RanP_ the 
odd sector. The following lemma is an easy application of definitions, and we give 
it without proof: 

Lemma 3.7. We have that tc{P±f,P^g) = for all f,g E n^{R). Then the 
domain Dc of He is left invariant by P± ; moreover, He commutes with P± . 

□ 

A standard consequence of the elliptic regularity (see also (|3.3U|) ') is the fol- 
lowing lemma, given again without proof: 

Lemma 3.8. The eigenvectors of He belong to C°°(M\{0}). 

□ 

A less obvious result is the following characterization of Dc'- 
Lemma 3.9. Every ^ in Dc obeys the following boundary condition at 0: 



lim 



21n(r)?A(0) - 21n(2£)V'(0) 



= 0. (3.32) 



18 



Proof. Now remember that tp' is continuous outside the origin if -0 G Dc, so ip'{zke) 
makes sense. Moreover, for every € 7i^(R) we have: 



- / (/.'(x)0'(x)da: + 21n(r)<^(O)0(O)~Co(0,0) = (0,i?cV'>- (3.33) 



i;'{x) = ^'(1) + / r{y)dy = ^'(1) - 2 / + U{y)) dy 



We can write for a; > (see (|3.30|l ') 

hi = w/n ^ - ■? 

/i V y 

and then for x close to we have: 

IV' (x)! < const + 2|lnx| HV'Hoo. (3.34) 

The same estimate is true for negative x near 0, and together with the estimate 
l|3.4|l it follows that ip' diverges at not faster than a logarithm. 
Now we can integrate by parts outside the origin and write: 



(j)' {x)^l^' {x)dx (3.35) 



£)i/;'(-e)- / (/)(x)i^"(x)dx - 0(e)V'(e) - / (l){x)ip" {x)dx 



4>' {x)'(Ij' {x)dx 



where the last integral will converge to zero with e. 
After a similar integration by parts we obtain: 



Co(</.>) = j H-2x) ■ {d^{H)){x)dx - j \n{2x) ■ {d,XH)){x)dx 
+ ln(2e) f 0(i)V'(e) + 0(-£)V'(-£)) + / A • l^iij{x)dx. (3.36) 



Following the reasoning in H3.18|l . one can prove that 



ln(2a:) • [d,{H)]{x)dx 0{ei\\ne\) Mn^Mln^ 



and thus 

a 



'o(0, V) ln(2£) Uie)^P{e) + 0(-e)V^(-£)) + / ^ • ^{x)ij;{x)dx. 

(3.37) 



19 



Putting \6.'S7\i and in H3.33|l . and using (|3.3Uj) . we eventually get: 



lim 



- HeWie) + 21n(r)0(O)V(O) 



ln(2£) (0(£)V(£) + 0(-e)V^(- 



0. 



(3.38) 



The last ingredient is the embedding H^(M) C C^/'^-^{R), and the estimate 
which allows us to simplify the limit: 



lim>(0) 

e— >0 



V''(-e)-^'(£) 



21n(r)V'(0) - 21n(2£)V'(0) 



for all (p G Ti.^{R). The lemma is proved. 



(3.39) 
□ 



3.2 The eigenvalues and eigenvectors of 



c 



In this subsection we give analytic expressions for eigenvalues and eigenvectors of 
He corresponding to the negative, discrete spectrum; much of the information we 
need about special functions is borrowed from (Lj . We want to have the same formal 
expression for our eigenvalue problem outside z = as in that paper, namely 



1 , a , 







Tr-2- Let us 



where '0 will be an eigenfunction with an associated eigenvalue E - 
now do this in a rigorous manner. 

We want to implement the change of variables x = ^az, a > 0, which 
amounts to defining a unitary operator on L^(R): 



2 V a 

Now assume is a normalized eigenvector for Hq satisfying 

Hc(t) = E(j), E < 0. 



'ijj{2z/a). 



(3.40) 



(3.41) 



Instead of solving the above equation, we will reformulate it in terms oi ip — Uw4>^ 
and (j) — U^ijj. To do that, we need to fulfill two conditions. The first one is: 



[UwHU^^^]{z)^^[-r{z)-^^P{z)]=Ei,{z), z^O, or 
r{z) - \m + ^Mz) = 0, z^Q,E^ 



(3.42) 



The second thing is to see what condition at z = should ip obey in order to be 
sure that U^ip is in the domain of He- If we replace tp in H3.32|l hy (p — U^ip, 
then we get the modified condition: 
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lim 



2 ^'(-le)-^y{le) 
a 2 



21n(r)i/'(0) - 21n(2e)V'(0) 



or 



lim 



+ a(lnr - ln(ae))?/'(0) 



= 0. 



(3.43) 



(3.44) 



Therefore we reduced the problem of finding the eigenfunctions and eigenvalues 
of He to solving the ordinary differential equations in 13.42|l . with the boundary 
condition given in (|3.44l) . We will see that solutions can be constructed only if 
a obeys some conditions. 

A priori, a can be any positive real number. First assume that a is not a 
positive integer. Then if we solve (|3.42|l for 2: > 0, we see that the only square 
integrable solution at +00 is the one given by a Whittaker function: 

M^„.i(z) = ze-^^C/(l - a,2,z), (3.45) 

where U is the confluent hypergeometric function or Kummer function, see |AS| . 
If a is a positive integer, the solution is obtained as the limit of 1 (z) when a 
tends to a positive integer TV and get: 



lim 1 (z) 

. AT o ^ ' 



(3.46) 



where L]^_^ is an associated Laguerre polynomial. 

We denote with r(z) and ip{z) — r'(z)/r(z) the usual gamma and digamma 
functions. We have the following first result: 

Proposition 3.10. (i). All negative eigenvalues of He are non-degenerate. The 
eigenfunctions of He are also eigenfunctions of P± . 

(ii) . There exists an infinite number of odd eigenfunctions (j)odd,k, k £ Z_(_, 
corresponding to every aG{l,2,...}. 

(iii) . There also exists an infinite number of even eigenfunctions 4'even,kj 
k G each corresponding to a certain € (fc — 1, fc) for every k G Z_|_. 

Proof, (i). Choose any eigenfunction (p of He corresponding to i? < 0. Make the 
change ip — U^4>^ and then look at the associated differential equation: 



V"(z)-i^(z) + ^V(^)-0, 



lim 



Q;(lnr — ln(Q!£))i/'(0) 



= 0. 



(3.47) 



(3.48) 



First assume a = N £ Z+. The theory of ordinary differential equations insure 
the existence of two constants Ci and C2 such that 

1 



) = C,e'-^^z-L\,_,{z) 



N 

C2e5^z-^L]v_i(-z), 



z > 0, 
z < 0. 



(3.49) 
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By inspection (and by continuity) we get that '0(0) = 0. If we put ip in 
and using the expUcit form of the Laguerre polynomials, we get that the boundary 
condition is fulfilled only if Ci — C2 which amounts to ip{z) = — z), i.e. there 
is one and only one solution which is also odd. 

Now assume a ^ Z+ . The theory of ordinary differential equations insure the 
existence of two constants C3 and C4 such that 

V'(z) = C3l^„,i(z), z>0, 

ipiz) ^ C'iW^ ii-z), z<0. (3.50) 

By inspection we see that 1 (0) = p^^^^-^ 7^ 0, hence by continuity at zero we 
must have C3 — C4, hence ip{z) = ip{—z) and the eigenfunction must be even. 

(ii) . The proof is already contained in (i), since the boundary condition is 
trivially fulfilled for odd functions. There is exactly one eigenfunction, an odd one, 
corresponding to every a € . 

(iii) . We saw in (i) that if there are eigenfunctions corresponding to a ^ 
then they must be even. In order to get all possible a's which are compatible with 
the boundary condition (|3.48() . we compute (note that ip' is odd): 



lim 



-W' 1 (e) + Q;(ln r - In(ae)) 



r(l - a) 



(3.51) 



and using the explicit expression of these special functions we obtain the condition 
on a: 

f{a,r) :='(/;(l-a) + 27H lna + lnr = 0. (3.52) 

2a 

where ip here means the digamma function and 7 is Eulcr's constant. Since the 
digamma function is strictly increasing from —00 to +00 on each interval of the 
form (— m, —m + 1), m g one can easily see that f{-,r) is strictly decreasing 
from +00 to —00 when a varies in an interval of the form (fc — l,k) for every 
k £ Z+. Therefore we have a unique solution ak G (A; — l,fc) of the equation 
f{ak,r) for every k £ Z^. The proposition is proved. 

□ 



The previous proposition stated that only the eigenvalues from the even sector 
can vary with r. Let us now further investigate this dependence. 

Corollary 3.11. (i). The excited states with even parity tend to those with odd 
parity when r is small. More precisely, for k > 2, we have that lim^^o olu — k~ 1; 
(ii). Fork = 1, we have the following asymptotic behavior of the ground state: 

«iW = -^l^{l + o,.(l)}, £;i(r) = -2[ln(r)]2{l + o,(l)}. (3.53) 
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Proof, (i). The limit follows easily from the properties of the digamma function, 
(ii). We apply the implicit function theorem. Define the function 



^ l + 2a[27 + V(l-a)-ln(a)] ^' 

for {a, y) in a small disk around the origin in M? . This function is near (0, 0), 
(9ctF)(0, 0) = 2, and F{Q, 0) = 0. Then for every y > small enough there exists 
a{y) > such that F{a{y),y) — 0. Now put y = — l/ln(r) and we are done 
because (13.52(1 is also satisfied with this a. □ 



3.3 Approximation of if^^ by He 

We will now show that the negative spectrum of H^g converges in a certain sense 
to the one of He- This is made precise in the next proposition, but before we need 
a definition. For a given subset S of M, and for any e > we define 

:= U B,ix). (3.54) 

If 5' is a discrete, finite set, then is a finite union of intervals of length 2e, 
centered at the points of S. 

Proposition 3.12. The following three statements hold true: 

(i) . Fix a < 0, and denote by A :— a{Hc) H (— cx3,a] and B := a[H^g) D 
(—CO, a]. Then for every e > 0, there exists r^ > such that for every r < r^ we 
have 

AdB^, Be A^. (3.55) 

(ii) . The ground-state of H^^ is non-degenerate, has even parity, and diverges 
to —oo when r 0. Moreover: 

lim I Ma{Hc) - Ma{H^ff)\ = 0. (3.56) 

(iii) . Fix a compact interval [a, b] C (— oo, 0) and suppose that He has exactly 
one eigenvalue of a given parity Ec in [a,b], for all r < r^. Then if r is small 
enough, H^^ has exactly one eigenvalue of the same parity E^s in this interval and 

\im\E,s~Ec\^0. 

Proof. Let us introduce the resolvents Rcs{z) = {H^ff — z)^^ for all z E p{H^ff) 
and Rc{z) — {He — z)^^ for all z G p{Hc). The first ingredient in the proof is 
contained by the following lemma: 
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Lemma 3.13. There exists a constant K > 1 sufficiently large, and ro small 
enough, such that for every r < r^ we have that the form defined on x 
(see also (|X77jl ) 

Vcif.g) tc {[vll2 + XrV^f, [pII2 + \r]-"^g) 
+ A,(/,b2/2 + A,,]-i.g>-(/,.g) 

" = "[p2/2 + K]-^'^{Hc + K}[pI/'2 + A.]-^/2 _ A, := K\n\r). 

(3.57) 

generates a hounded operator on L^(K) denoted in the same way. Moreover, supq^^^^^^ I I^cI I < 
1/2. And we have 

{He + AJ-i = [pI/2 + A,]-i/2{Id + Vc}-\pI/2 + K]-^l\ (3.58) 
Proof. The key estimate is contained in 

\\{pll2 + A)-i/2|U.^i» < const (3.59) 

obtained with an argument as in (|3.15|) . Then if we have | ln(r)|/v^ small enough, 
then the "delta function" part of tc will be small uniformly in r < tq. Using also 
l|3.19|l . and the definition (|3.27|l . then one can show that (|3.57|l is a bounded 
sesquilinear form on L^(R), with a norm which can be made arbitrarily small if K 
is chosen large enough. Now the equality (|3.58|l is easy, and note that this is also 
compatible with (|3.53|) . □ 

Introduce the notation: 

yeff {pI/2 + Kr''^VM/2 + Xrr'/\ (3.60) 

The second ingredient in the proof of the above proposition is the following esti- 
mate, which is an easy consequence of Proposition 13. II 

||^c-V;ff|| = 0(r4/9), r<ro. (3.61) 

We also have that ||T4ff|| < 2/3 if ro is small enough, uniformly in r < ro, and 
then 

{Hls + \r}-^ - [pl/2 + \r]-^/^{ld + V,s}-\pl/2 + \r]-^l\ (3.62) 

It is clear that a similar identity would hold for any other A > A^, and this already 
tells us that the spectrum of H^^ is contained in an interval of the type (— A^, oo), 
thus justifying the discussion after 1)2.14(1 . 

From (|3.62l) . 13.61|l and (|3.58|) . we get the crucial estimate: 

4/9 

||i?eff(-A^)-i?c(-Ar)||< const r<ro. (3.63) 
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This estimate allows us to prove (i). Introduce the notation 

dc{z) := dist(z, a{Hc)). (3.64) 
Choose z e p{Hc) (thus dc{z) > 0). From the identity: 

{He - z)Rc{-Xr) = Id - (z + Xr)Rci-\r) (3.65) 
we get that the right hand side is invertible and: 

{Id - (z + X)Rc{~Xr)r^ - {He + Xr)Rc{z) = Id + (z + Xr)Rc{z). (3.66) 
The first equality implies that 

Rc{z) = Rc{~K) {Id - (z + Xr)Rc{-K)y^ , (3.67) 
while the second one gives the norm estimate: 

{Id - (z + Xr)Rc{~Xr)y^ II - + (1^1 + Xr)/dc{z)]. (3.68) 

Note the important fact that H3.67|l is just another form of the resolvent identity, 
valid for any self-adjoint operator. If we could replace Rc{—Xr) by Rcs{—Xr), then 
the right hand side would immediately imply that z G p{Has). 

We can restrict ourselves to those z's which obey |z| < A^. Then using H3.63|l 
and (|3.68|) . we get that for r < rp and dc{z) > r^/^, the operator 

Id - {Z + Xr)Rcs{-Xr) 

is invertible and we get the estimate: 

{Id - (z + A,)i?eff(-A,)}"' < Xr [1 - const A, r^/ydciz)]-\ (3.69) 

dc{z) 

Therefore we have proved that for every z which obeys |z| < A^ and dciz) > r^/'^, 
the operator 

i?eff(-A,) {Id - (Z + Xr)Ros{~Xr)r^ (3.70) 

exists and defines Rcs{z)- It means that the spectrum of H^g is "close" to that of 
He, and the distance between them is going to zero at least like r^/'^. 

Let us now prove (ii). We know that the ground-state of He diverges like 
— In^(r) for small r, and it is isolated from the rest of the spectrum. Choose a 
circular contour T of radius 1 around this ground-state. It means that dciz) — 1 
for z e r. Then jXTUIl . and imply the estimate 

sup ||i?eff(z) - i?c(2)|| < const -r^/^, r < tq. (3.71) 
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Now we can employ the regular perturbation theory, see by using Riesz pro- 
jections defined as complex integrals of the resolvents on contours like F. Then the 
estimate (|3.56|l is straightforward. 

Finally, let us prove (iii). We know that for small r, the excited states of He 
tend to cluster in pairs. The eigenvalues from the odd parity sector are independent 
of r, while those from the even parity sector will converge from above to the odd 
ones (see Proposition 13 . 10|l . Consider such a pair of eigenvalues, which will always 
remain separated from the rest of the spectrum if r < ri and ri is small enough. 
Then we can find a contour F which contains them and 'miz^v dc{z) is bounded 
from below uniformly in r < ri . Then we can again write an estimate like (I3.71|) , 
and then apply the regular perturbation theory. The proof is over. 

□ 



4 Reduction of to 

We are now ready to go back to (|2.24l) , and argue why only the diagonal entries of 
the infinite operator-valued matrix {Hm.n}m,n<£Z are important for the low lying 
spectrum of H^. 

Let us formally write as: 

= i?diag + Kjffdiag, 

where Hdiag — ®„gz(-^cff + ip-), and yoffdiag contains all the non-diagonal entries 
of the form „, m 7^ n, (see (I2.25|) '). and zero on the diagonal. We will prove 
in this section that T4ffdiag is relatively form bounded with respect to -ffdiag, and 
moreover, it is a "small" perturbation when r is small. 

The main result is very similai^to Proposition 13 . 1 31 where we only have to 
change He by i/diag and H^g by if. Parity here only refers to the x variable. 
Therefore we will start comparing the two operators. Before that, let us note that 
the negative spectrum of i?diag is given by H^g if r is small enough. 



4.1 Kflfdiag is ifdiag-form bounded 

Let Xr = K\n^{r) with K large enough and r < tq, as in the previous section. We 
know that — G p(^^diag), and denote by i?diag(— Ar) the resolvent (i?diag + Ar)~^- 
Then the main technical result of this subsection will be the following estimate: 
there exists 6 > and ro{6) such that 



RdiLi-^r)Vomi.sR'Ji!j-Xr) ^0{r'X-'/^), r < tq. (4.1) 

aiag V / & aiag V / ^ (Z-L^ (""^ ' ' \ f v / 



Let us first notice that we can replace i^^iag ^ simpler operator, namely ®„gz(ep^" 
1 -|- ^^)~^^'^, where e is a small enough positive number. Indeed, we can write 

(/, [H^g + Xr + nV(2r2)]/) > (/, [epl + A./2 + nV(2r2)]/), 
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< 1. 



(4.2) 



where we used that for e smah enough we can show that: 

(1/2 - e)p2 _ 1^;^. + A,/2 > 0, r<ro. 

This means that 

[epl + \r/2 + nV(2r2)]i/2[iy;^ + a. + ny{2r')]-^/' 

Define the bounded operators in i^(R) (see (|2.25ll ): 

:= [^pI + A./2 + my{2r')]-'/'v:,Jepl + A./2 + {2r')]-'/' , m + n, 
V;;,^:=0, meZ. (4.3) 

Then (|4.1|) would be imphed by the fohowing, stronger estimate: 



I ^ m.n ( 



ro- 



(4.4) 



B„(i2(Z:L2( 

By an easy apphcation of the Schur-Holmgren lemma, one can prove the estimate: 



I * m,n f 



< sup \\V^,n\\B^(L^ 



(4.5) 



We now concentrate on the norms \ \V^ nl |boo(L2(r)) and study their behavior 
in r, m, and n. Remember that only the case m ^ n is of interest, since the 
diagonal terms are zero. 

Before anything else, let us do a unitary rescaling of L^(R) by {U f){x) := 
r^^^f{rx) and {U*f){x) := r^^/"^ f{x/r). Then due to various homogeneity prop- 
erties we get: 

t^C^nt^* = r ■ [epl + r'Xr/2 + mV2]-i/V^,„[ep2 ^ ^ly^j^ ^ /2]-'/\ (4.6) 
We first give an important estimate for „, stated in the next lemma: 

Lemma 4.1. Let < a < 1 and \m — n\ > 1. Fix any < e < 1. Then there 
exists a constant C — C(a,e) such that we have the following estimate: 



\V'\{x) < C 



1 



1 



1 



\n — mr \x 



and 



l^m,nl(a;) < const - 



^3e 



l^l > r 



(4.7) 



(4.8) 
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Proof. Due to symmetry properties we can write 



Integrating by parts we get: 

.ri , . 1 r sin[(» - m)y] sin(y) 

"""-("^ = M^T^^ i-. [.^ + 4sin^(,/2)]3/2'^- 

This equality immediately proves H4.8|l . So we now focus on \x\ < r~'^. We can 
split the integral in two: one in which \y\ > 7r/2, and where the integrand has no 
singularities when x is small, and the second where \y\ < 7r/2. In that region we 
can use the same idea as in Lemma 1301 of replacing sin^(y/2) by y^. We hence get: 

\n-m\y J_^/2 [x^+y^]y^ J 
Now we employ the inequalities (here < a < 1 is arbitrary): 

|sin[(n-m)y]| < |n - m|i^"|y|i'", |sin(y)| < \yl 
then we make the change of variables s — y/\x\ and write: 

/o [1 +^3/2 

Thus the lemma is proved. 



, const / „|n-TO|i-" f°° , , 

VlJix)<j-—^il + 2 ^ / -——ds). (4.12) 



□ 



Now let us go back to (14.61) . and estimate the various norms. If we write 
Kl,n = VA,nX{\ ■ I < + Kl,«x(l • I > r~'), then we have two different types of 
estimates. When we keep „x(| • | > r"*^), which is bounded, then for the two 
resolvents we can use the usual Boo{L^) norm, which together with 1)4. 8|l gives a 
contribution: 

const r^*^ , / , , „x 

n^m. (4.13) 



y/r^Xr + ^r^Ar + \m — n\ 

When we keep „x(| • | < r''^), the estimate from (I4.7|l gives us that 
[iKn.nlxd ■ I < 1)]^^^ is an function, hence the operator 
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is Hilbert-Schmidt. Thus we have a product of two Hilbert-Schmidt operators, and 
we can give an upper bound for the Boo norm of their product of the form: 

const(Q!) r"^ 

^ ' n^m. (4-14) 



Therefore we obtained an upper bound for the norm of the operator in l|4.6|l 
of the form: 



m,n I 



r ■ const • r~'^ 
(r2A,. + n2)i/4 [r'^X^ + ^2)1/4 |^ _ „|a ■ 



m^n. (4-15) 



Remember that one is interested in the right hand side of 1)4. 5|l . Now choose 
1/2 < a < 1. We have to investigate several cases: 

1. When m — and |n| > 1. Then the first term in H4.15|l wiU behave hke 

The second term will behave like r^/2-«^,^ ^^^|n|^^/2-a Both contributions 
are summable with respect to n. Note that if e is small enough, both expo- 
nents of r are positive. Denote by S the smaller one. 

2. Fix m ^ 0, and consider all n ^ m. When n = 0, we get similar terms as 
above. If n ^ 0, then we remain with the problem of summing up something 
hke ^ 

sup |m|~"^/^ , ,,,„, ;— , 1/2 < a. 

We can either use Holder's inequality, or we can split the above sum in the 
following way: 

y I 

|n|i/2|n-TO|" 

n^O,n^m III I 

E + E 

I n^O,n^?-n,|n| < |n— m| n^0,n^7n, |n| > |n — m| y 

< I I I L /9 + 1 TTT^ ) < const(a). (4.16) 

We therefore consider H4.1|l as proved. 
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4.2 Comparison between and -ffdiag 
If r is small enough, we have the identity: 

(ff- + A.)-i = i?^Lg(-Ar) {id + i?d£(-A.)Kffdiagi?^£(-A.)}"' i?d(4(-A.). 
Moreover, this implies: 

IKiF + A,)-l - i?diag(-A„)|| < £^||i?V2^(-A,,)Kffdiagi?d£(-A.)|| 

< const — ^ . (4.17) 

This is the same as what we had in (I3.63f) . but with A^ instead of r*/^. 

Therefore we can repeat the arguments of Proposition 13.121 and prove a similar 
kind of spectrum stability for and -ffdiag- 



5 The main theorem and some conclusions 

We now try state a concentrated main result of our paper. Let us first go back to the 
very first Hamiltonian which was declared to be relevant for the exciton problem. 
This is , written in (|2.1|) . Because of the Coulomb singularity, the best way to 
look at the spectral problem is to consider its form tn, given by (|2.4|l . H2.5|l and 
H2.6|l . We then managed to separate the mass center motion in the longitudinal 
direction, and we got a simpler form ty^ in H2.8II . The center of the mass cannot 
be separated in the transverse direction because of the cylindrical geometry, but 
at least we can write th as a direct sum of ^^.^^ t/i^ . A crucial observation has 
been stated in (|2.14|l . which says that only th^ is responsible for the lowest lying 
spectrum of the original form. 

This gave us the possibility of renaming th„ with tu in H2.15|l . and declare 
it as the central object of study. Then in Proposition 12 . II we constructed its asso- 
ciated self-adjoint operator , where we had to take care of the Coulomb-type 
singularity in two dimensions. 

Then after a unitary transformation induced by the discrete Fourier transform 
with respect to the y variable, we can see as an infinite operator valued matrix 
acting on the Hilbert space ^^(Z;L^(R)). We then decomposed as the sum of 
a diagonal operator H^is^g and an off-diagonal part V^ffdiag • Eventually we proved 
in Section 0] that the low lying spectrum of is only slightly influenced by the 
off-diagonal part for small r, and therefore the relevant object remains i?diag- 

But this diagonal part has the nice feature that each of its entry is of the 
form H^g + ^,neZ, where ff^'^fj is given in fTI^ and (|T77I) . Then in SectionEl 
more precisely in Proposition l3.12l we prove that the low lying spectrum of H^g is 
well approximated by the spectrum of a solvable operator. He, which we discussed 
in Proposition 13 . 1 01 
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We are now ready to collect all these results in the main theorem of our 
paper: 

Theorem 5.1. The following three statements hold true: 

(i) . Fix a < 0, and denote by A :— a{Hc) H (—00, a] and B :— a{H^) D 
(—00, a]. With the definition introduced in (|3.54|l . we have that for every e > 0, 
there exists > such that for every r < r^ we have 

AcB,, Be A,. (5.1) 

(ii) . The ground-state of is non- degenerate, and diverges to —00 when r — > 
0. The corresponding eigenfunction has even parity with respect to both variables. 
Moreover: 

lim I inf o-(iJc') -mfo-(ir'-) I = 0. (5.2) 

(iii) . Fix a compact interval [a, b] C (—00, 0) and suppose that He has exactly 
one eigenvalue Ec in [a,h], of parity p — ±, for all r < tq. Then if r is small 
enough, has exactly one eigenvalue E in this interval and 

lim 1^ - Ec\ = 0. 

Moreover, the corresponding eigenfunction has parity p with respect to x. 

Another important aspect of this problem is to determine how fast these 
limits are assumed. We have not touched this issue here, but we will study the 
numerical and physical implications of our results in a consequent paper. 
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